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ABSTRACT 

High resolution N-body simulations of cold dark matter (CDM) models predict 
that galaxies and clusters have cuspy halos with excessive substructure. Observations 
reveal smooth halos with central density cores. One possible resolution of this conflict 
is that the dark matter is warm (WDM); this will suppress the power spectrum on 
small scales. The Lyman-alpha forest is a powerful probe of the linear power spectrum 
on these scales. We use collisionless N-body simulations to follow the evolution of 
structure in WDM models, and analyze artificial Lyman-alpha forest spectra extracted 
from them. By requiring that there is enough small scale power in the linear power 
spectrum to reproduce the observed properties of the Lyman-alpha forest in quasar 
spectra, we derive a lower limit to the mass of the WDM particle of 750 eV. This limit 
is robust to reasonable uncertainties in our assumption about the temperature of the 
mean density gas (To) at z = 3. We argue that any model that suppresses the CDM 
linear theory power spectrum more severely than a 750 eV WDM particle cannot 
produce the Lyman-alpha forest. 

Subject headings: cosmology:theory, dark matter, methods:numerical, quasars: 
absorption lines 

1. Introduction 

Cold dark matter (CDM) models of structure formation have been remarkably successful in 
explaining a wide range of observations at both low and high redshifts. However, high resolution 
N-body simulations of the CDM model have revealed two potential conflicts with observations. 
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First, the inner mass density profiles of simulated halos are cuspy, as opposed to the shallow 
profiles inferred from the rotation curves of dwarfs and low surface brightness galaxies ( pVloorej 
1994| |de Blok fc McGaugh 1997| ; [Hernandez fc Gilmore 1998| ; but see van den Bosch et al. 1999), 



and the observations of soft cores in galaxy clusters by gravitational lensing ( Tyson, Kochanski,| 
fe DellAntonio 1998Q . Second, group-sized halos in simulations contain a large number of low 
mass subhalos, greatly in excess of the observed number of satellite galaxies in the Local Group 
QKlypin et al. 1999| ; |M oore et al. 1999| ), A variety of mechanisms have been proposed recently 
to reduce the clustering of mass on small scales, while simultaneously retaining the large scale 
successes of the CDM model. These mechanisms include breaking the scale invariance of the power 
spectrum of primordial density fluctuations ( Kamionkowski Sz Liddle 1999| ), postulating different 



collapse histories for dark matter and baryons (Navarro, Eke, Sz Frenk 1996; Bullock, Kravtsov, 



fc Weinberg 2000| ; Binney, Gerhard, &: Silk 20001) , and modifying the nature of dark matter. The 
last option includes dark matter models with a variety of properties — self-interacting ( |S per gel 



fc Steinhardt 2000| ), warm ( [Sommer-Larsen &: Dolgov 1999| ) repulsive ( Poodman 2000| ), fluid 
(Peebles 2000| ), and fuzzy ( [Hu, Barkana, & Gruzinov 200q) . 

In the warm dark matter (WDM) model, the linear power spectrum is exponentially 
damped on scales smaller than the free-streaming length of the warm particle, 
R f = 0.2(O w /i 2 ) 1/3 (m w /keV)- 4 / 3 Mpc, relative to the pure CDM model ( [Bardeen et al 



1986 ). Here, f2\y is the ratio of WDM density to the critical density, mw is the mass of the 
WDM particle, and h = -ffo/100 km s _1 Mpc -1 is the Hubble parameter. Non-linear gravitational 
evolution transfers power from large scales to small scales, so the clustering of the highly nonlinear 
mass distribution is insensitive to the shape of the linear power spectrum below the non-linear scale 
([Little, Weinberg, fc Park 1991fc |Bagla fc Padmanabhan 19971 |White fc Croft 200q ). Conversely, 
Ly-a absorbers seen in the spectra of high redshift (z ~ 3) quasars arise from mass density 
fluctuations in the quasi-linear regime ( |Bi 1993 ; Hernquist et al. 1996; Hui & Gnedin 1997), so 
their properties remain sensitive to the linear power spectrum. In this Letter, we set a lower limit 
on the mass of WDM particle by requiring that there be enough small scale power in the initial 
conditions to reproduce the observed properties of the Ly-a forest at z = 3. 



2. Simulations 

We study the Ly-a forest in both CDM and WDM models with VL m = 0.4, Q A = 0.6, 
h = 0.65, n b h 2 = 0.02, and &8m = 0.95, where fi m , f^A, and £lb are the contributions from total 
mass, vacuum energy, and baryons to the total energy density of the universe, and as m is the rms 
density fluctuation in 8/i -1 Mpc spheres, chosen here to reproduce the observed cluster abundance 



( |White, Efstathiou, fe Frenk 1993| ; |Eke, Cole, fc Frenk 1996| ). Table 1 lists the parameters of all 



models investigated. Our WDM models have mw = 200, 500, 750, and 1000 eV (corresponding to 
Rf = 0.95,0.28,0.16, and O.llMpc), spanning the range of WDM masses required to match the 
phase space density of dark matter cores in disk galaxies and dwarf spheroidal galaxies ( [Hogan fc 
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Dalcanton 2000 and references therein; |Sellwood 2000 ) We also analyze a Broken Scale Invariance 



model (BSI, [Kamionkowski fc Liddle 1999| ), using an analytic fit to its power spectrum from White 
&i Croft (2000), with a filtering scale ko = 2h Mpc -1 required to fit the observed number density 
of low mass satellite galaxies in the Local Group. 

We calculate the linear power spectrum of the mass density field for all the models using 
the full Boltzmann transport code of Ma & Bertschinger (1995). We assume the usual three 
massless neutrino species with a present-day temperature Tq jV = 1.947 K, and treat the WDM 
component as a fourth (massive) species with Tq^/ = Tq u (^V/i 2 ) 1/3 (93eV/mvi/) 1/3 . Fifty 
Legendre moments are used to follow the evolution of the WDM phase space. We compared our 
WDM power spectra with the fitting function of Bardeen et al. (1986) and find that their formula 
provides a reasonable approximation for k < 5h Mpc -1 [if a baryon-corrected shape parameter 
T = Q m h exp(— (1 + fi^, )) is used], but it overestimates the power spectrum by up to an order 
of magnitude at higher k. 

We employ a particle-mesh (PM) N-body code that is described in detail in Steed et 
al. (in preparation). This code computes forces using a staggered mesh (Melott 1986; Park 



1990D , and integrates the equations of motion using the leapfrog scheme described in Quinn 
et al. (1997). Our periodic simulation volume has L = 25/i _1 Mpc, with N p = 256 3 particles 
and an N m = 512 3 force mesh. We assign initial displacements and velocities to the particles 
using the Zeldovich approximation. In the WDM models, we also add a randomly oriented 
streaming velocity component that is drawn from a Fermi-Dirac distribution with an rms velocity, 
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v ims = 48.7 (gAg) 3 (io5^v)~ (ft) 3 km s_1 - We evolve from redshift z = 49 -» 3 in 24 

equal steps of the expansion scale factor. 

Figure 1 shows the linear and non-linear power spectrum of the mass density field at z = 3 
in different models. The differences between different models are significantly smaller in the 
non-linear power spectra, compared to the differences in their linear theory power spectra, because 
non-linear gravitational evolution has regenerated power on small scales. Nevertheless, power 
transfer is not entirely successful in erasing the differences between the power spectra of different 
models. Thus, at z = 3, the WDM200 model has almost an order of magnitude less power on 
small scales compared to the CDM model, while the models WDM500, WDM750, and BSIK02 
are all deficient by ~ 0.3 dex. We also ran these models with a 50/i _1 Mpc box using the same N p 
and N m , and found similar results on overlapping length scales. 

The model CDM200 has the same linear power spectrum as WDM200, but v Tms = 0. This 
model gives virtually identical results to the WDM200 model, showing that the velocity dispersion 
of WDM particles has a negligible effect on the Ly-a forest at z = 3. Hence, we do not show the 
results for this model in the figures. 
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Fig. 1. — Power spectra of the mass density fields in different models, (a) Linear mass power 
spectra at z = 3 from the Boltzmann code, (b) Non-linear power spectra at z = 3 from simulation 
output, computed by cloud-in-cell binning onto a 256 3 grid, taking its Fourier transform, and 
squaring the amplitudes of its Fourier components. Note that WDM750T025K has the same mass 
power spectrum as WDM750. 
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3. Properties of the Ly-a forest 

The Ly-a forest arises from a continuous, fluctuating, low-density [— 1 <^ log(p g /p g ) <J 1] 
gaseous medium whose temperature is determined by the balance between photoionization heating 
by the UV background and adiabatic cooling due to Hubble expansion (see e.g., pi 1993j ; |C cn 



|et al. 1994 ; Zhang, Anninos, fc Norman 1995 ; Hernquist et al. 1996| ; [Hui Gnedin 1997 ). 



Photoionization equilibrium leads to a power-law relation between temperature and gas density 
of the form T = T (p g /p g ) a , with 4000 K < T < 10000 K, and 0.3 < a < 0.6 (|Hui fc Gnedin 



1997 ). The optical depth to Ly-a photons (TL yQ ,) is proportional to H I gas density, which in turn 
is proportional to T~ ' 7 for gas in photoionization equilibrium near 10 4 K. Pressure gradients 



in this low-density, fairly cool (T ~ 10 4 K) gas are small, so gas traces dark matter ( Proft et 



al. 1998| ). Hence, the gas density-temperature relation can be translated into a relation between 



mass density (p m ) and TL yQ (the "Fluctuating Gunn-Peterson Approximation", Rauch et al. 1997] ; 



|Croft et al. 1998] ; |Weinberg, Katz, fc Hernquist 199S| ), Thus, the transmitted flux of Ly-a photons 
[F = exp(— TL ya )] in a quasar spectrum is a map of the mass density fluctuations along the line of 
sight. 

We extract artificial Ly-a absorption spectra along 400 random lines of sight through 
each simulation, using the TIPSY package ( [Katz h Quinn 1995| ; see [http:// wwwH 



hpcc.astro.washington.edu / tools /TIPSY ) . TIPSY calculates the local mass density (p m ) at 
the position of each dark matter particle using a cubic spline smoothing kernel (Hernquist & Katz 
1989) enclosing 32 neighbors, and assigns it a temperature T = To(p m /p m )°' e . Spectra are then 
computed from this particle distribution using the algorithm described by Hernquist et al. (1996). 
We sample the spectra with resolution Av = 6.4 km s , roughly corresponding to that of the 
Keck HIRES spectra to which we will compare our model predictions. We "fit a continuum" to 
each spectral segment, 3200 km s _1 in length, by rescaling all flux values so that the highest pixel 



has transmission of unity ( [Dave et al. 1999 ). Fitting a continuum makes low-power models appear 



more like high power models, so it is an important step for conservatively discriminating against 
models with suppressed power. 

The optical depth to Ly-a photons depends on a number of uncertain parameters including 
rife, h, Tq, and the intensity of the ionizing background. However, the cumulative effect of all these 
parameters is to scale the normalization factor (A) relating TL yQ to p m /~p m , as long as collisional 
ionization is unimportant ( Weinberg et al. 1998[ ). We fix A in each model so that the mean 



transmitted flux at z = 3, F = (exp(— TL ya )), is 0.684, as measured by McDonald et al. (1999) 
from Keck HIRES spectra of eight quasars. 

Figure 2 shows artificial Ly-a absorption spectra at z = 3 extracted along four random lines 
of sight, in three different models: CDM, WDM500, and WDM750. The gas responsible for Ly-a 
absorption is in larger and more diffuse structures in WDM models compared to the CDM model. 
Hence, the WDM spectra are smoother and have less small scale power as compared to the CDM 
spectra. The fluctuations in the absorbing gas are also weaker in the WDM models, resulting in 
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absorption lines that are broader and shallower than corresponding lines in the CDM artificial 
spectra. 

Figure 3 shows the one dimensional power spectrum of the transmitted flux, Pp(k), computed 
using 400 artificial spectra for each model. The solid points show Pp(A;) measured by McDonald 
et al. (1999). The shape and amplitude of Pp(k) of the CDM model matches the observations 
very well, in the entire range —2.5 <^ log[/c(km -1 s)] <J —1. However, the WDM models predict a 
wrong shape of Pp(fc). Thus, although the Pp(^) of WDM models match the observations on large 
scales, —2.5 <^ log[A;(km _1 s)] <^ —2, they have too little power on smaller scales. This confirms the 
visual impression from the spectra in Fig. 2, and is also evident in the flux decrement correlation 
function (not shown). A higher value of To further suppresses Pp{k) on smaller scales, worsening 
the discrepancy with observations. Note that lowering O m to make up the power suppressed by 
WDM is not feasible; a lower Q m raises the power spectrum on all scales at z = 3, and thus would 
be inconsistent with the observed Pp(k) on large scales. 



Figure 4 shows the probability distribution function (PDF) of the transmitted flux (Jenkins 



O striker 1991; Miralda-Escude et al. 1996; Rauch et al. 1997[ ) , computed using 400 artificial 



spectra for each model. The solid points show the flux PDF measured by McDonald et al. (1999). 
The flux PDF of the CDM model matches observations remarkably well over the entire range of 
flux values. As the mass of WDM decreases, there is less small scale power in the mass distribution, 
the mass fluctuations are weaker, and the gas responsible for Ly-a absorption is in more diffuse 
structures. Hence there are fewer pixels where the transmission approaches unity, and more pixels 
at intermediate fluxes. In the flux range 0.6 < F < 0.8, the flux PDF for all WDM models with 
mw < 750 eV differs from observations by ^ 3a. This mismatch further confirms that the mass 
fluctuation amplitude of these WDM models is too low at small scales when the fluctuations 
are normalized to match Pp(fe) on large scales. Even the model WDM750T025K, with a higher 
temperature for the mean-density gas (To = 25, 000 K), does not match the observations as well as 
the CDM model. On the other hand, while WDM1000 does not match the observations very well, 
a higher value of Tq can bring its flux PDF to better agreement with observations. Since recent 



determinations of To prefer values around 20,000 K (e.g., pchaye et al. 1999|) , we conclude that 
the observed flux PDF cannot be reproduced by WDM models with m\y <^ 750 eV. The BSIK02 
model is also inconsistent with observations, being similar to WDM500. 



4. Conclusions 

For a given cosmology, the properties of the Ly-a forest depend on three factors: (a) the 
linear theory power spectrum, taken from a Boltzmann code; (b) power transfer from large to 
small scales during non-linear gravitational evolution, followed using N-body simulations; and 
(c) the equation of state of the intergalactic medium, governed by photoionization equilibrium. 
Here, we have used the observed Ly-a forest at z = 3 to constrain the shape of the linear power 
spectrum, and consequently m\y- Our principal conclusions are as follows: 
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Fig. 2. — Examples of artificial spectra at z = 3 along four random lines of sight in the simulations. 
All simulations have identical phases for the Fourier components of the initial density fields. In each 
panel, black line shows results for CDM, blue line shows WDM500, and red line shows WDM750. 
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Fig. 3. — Power spectrum of the transmitted flux at z = 3 from 400 artificial spectra for each 
model. Solid points with error bars show the flux power spectrum at z ~ 3 measured by McDonald 
et al. (1999) from eight Keck HIRES quasar spectra. 
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Fig. 4. — Probability distribution functions (PDF) of the transmitted flux at z = 3, computed 
using 400 artificial spectra for each model. Solid points with error bars show the flux PDF at z = 3 
measured by McDonald et al. (1999) using Keck HIRES spectra of eight quasars. 
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(1) m\v ^ 750 eV is required to match Pp(fc) and flux PDF of artificial spectra extracted from 
simulations of WDM models with Keck HIRES observations of McDonald et al. (1999), for 
any reasonable Tq. This limit is similar to that obtained from the maximum observed phase 
space density of dwarf spheroidal galaxies (Dalcanton & Hogan 2000). High resolution N-body 
simulations of WDM models with m\y ~ 1 keV predict the observed circular velocity function of 
satellite galaxies in the Local Group, but produce cuspy halos ( polin, Avila-Reese, fc Valenzuela 
120001) . 

(2) Ly-a forest properties depend on the shape of the linear power spectrum, and are insensitive 
to streaming velocities of WDM particles. Even for a 200 eV WDM particle, i> rms = 1.5 km s _1 at 
z = 3, an order of magnitude smaller than the thermal velocity of H I gas at T ~ 10 4 K. Hence, 
any mechanism that suppresses the CDM linear power spectrum more severely than a 750 eV 
WDM particle will be inconsistent with Ly-a forest observations. 

(3) The BSI model proposed by Kamionkowski & Liddle (1999), with ko = 2h Mpc -1 chosen to 
match the abundance of dwarf galaxies, does not reproduce the observed properties of the Ly-a 
forest at z = 3. We note that this is the extreme BSI model considered by White & Croft (2000), 
hence our conclusion is not in conflict with theirs. 
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Table 1: Model parameters. 



Model 


Power spectrum 


mw 


"rms(^i — 49) 


T (z = 3) 






(eV) 


( km s _1 ) 


(K) 


CDM 


CDM 







6000 


WDM1000 


WDM 


1000 


2.3 


6000 


WDM750 


WDM 


750 


3.3 


6000 


WDM500 


WDM 


500 


5.7 


6000 


WDM200 


WDM 


200 


19.2 


6000 


CDM200 


WDM* 


200 





6000 


WDM750T025K 


WDM 


750 


3.3 


25000 


BSIK02 


CDM^Slt 







6000 



"This model has the same linear power spectrum as WDM200. 

*We use an analytic fit to the power spectrum of the Broken Scale Invariance model with ko = 2h Mpc -1 , given by 
White & Croft 2000. 



